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Abstract 

We discuss the Penner-type matrix model which has been proposed to explain the AGT 
relation between the 2-dimensional Liouville theory and 4-dimensional M = 2 superconformal 
gauge theories. In our previous communication we have obtained the spectral curve of the ma- 
trix model and showed that it agrees with that derived from M-theory. We have also discussed 
the decoupling limit of massive flavors and proposed new matrix models which describe Seiberg- 
Witten theory with flavors Nj = 2, 3. In this article we explicitly evaluate the free energy of 
these matrix models and show that they in fact reproduce the amplitudes of Seiberg-Witten 
theory. 



1 Introduction 



Recently a very interesting relation between the Nekrasov partition function of M = 2 conformal 
invariant SU(2) gauge theory and the conformal block of the Liouville field theory was proposed 
[I]. It seems that this is the first example of a precise mathematical relationship between 
quantum field theories defined at different space-time dimensions. There have been various 
attempts at checking this AGT relation at lower instanton numbers by direct evaluation of 
Liouville correlation functions (2J |3l HI El [6]. There have also been attempts at proving the 
relation by comparing the recursion relation satisfied by the descendants of the conformal 
blocks and Nekrasov's partition function [TJ, [HI El EH] • 

On the other hand, a Penner type matrix model has been proposed to interpolate between 
the Liouville theory and gauge theory [XT] and provide an explanation for the AGT relation. 
In a previous communication [12] we have studied this matrix model and also proposed models 
for asymptotically free theories obtained by decoupling some of massive flavors. We have 
shown that the spectral curves of these matrix models reproduce those based on the M-theory 
construction and their free energies satisfy the scaling identities known in the SU(2) Seiberg- 
Witten theory. (See also [131 EH] f° r A r quiver matrix model). 

In this paper we would like to evaluate the free energies of these matrix models in the 
large N limit explicitly and show that they in fact exactly reproduce the amplitudes of SU (2) 
Seiberg-Witten theory. 

In section 2 we first describe the general properties of matrix models. In section 3 we 
compute the free energies: we integrate the Seiberg-Witten differential of the matrix model 
and evaluate the filling fraction in terms of the parameters of the spectral curve. We then 
invert this relation and derive the free energy. We present the computation for SU(2) gauge 
theory with two, three and four flavors and show that they all reproduce the amplitudes of 
Seiberg-Witten theory. In section 4 we discuss decoupling limits of some quiver gauge theories. 
Section 5 is devoted to conclusion and discussion. 

Note: in our convention the free energy of the matrix model F m is off by a factor 4 from 
that of gauge theory. Thus we will check the agreement 4F m = F gauge throughout this paper. 

2 SU(2) gauge theories and matrix models 

It has been proposed that the Nekrasov partition function for M = 2, SU(2) gauge theory 
with four flavors (summarized in appendix [Aj coincide with the four-point conformal block of 
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Liouville theory [TJ: 

ZZ {2) = Zcft = (^(00)^(1)^(^)^(0)) . (2.1) 

Here Vm is the vertex operator, Q = b + 1/b and the central charge of the Liouville theory is 
c= 1 + 6Q 2 . 

In order to relate the Liouville theory to matrix model, we consider the Dotsenko-Fateev 
integral representation of the four-point conformal block in terms of the free field <j)(z) [IB"]: 

Z DF = (|(/ dX i : ebHXl) : ) ^=0(00)^1(1)^,(9)^0(0)^ > (2-2) 

where the vertex operator VmX z i) i s given by : e mi ^ Zi ' : and we have introduced the iV-fold 
integration of screening operators. OPE of the scalar field is given by <fi(z)(j)(uj) ~ —2 log(z — oj). 
Momentum conservation condition relates the external momenta and the number of integrals as 
Y^=o ™i + + = Q. We redefine the momenta as m; = ^ for i = 0, 00 and rhi = ^ + ^ 
for i — 1, 2 pQ. Then the above condition becomes 

2 

^ irrii + tmpo + 2bg s N = 0. (2.3) 

As pointed out in [T6], [TT] and recently in [TTJ in the context of the AGT relation, the Dotsenko- 
Fateev representation may be identified as the /3-deformation of a one matrix integral 

Z DF = q^(l - q)^ (f[ I d\\ n(A, - \j)-^e^ W ^\ (2.4) 

\J=1"' / KJ 

In the case of b — i, integrations over {A/, I — 1, ■ ■ ■ , N} becomes an integral over a hermitian 
matrix M with eigenvalues {A/} and the action 

2 

W(M) = Y, m i M M - *)> z ° = 0> z ± = ^ = <?• (2.5) 

We identify the parameters mj with the mass parameters of the corresponding gauge theory. 
The identification of the parameter b with the Nekrasov's deformation parameters is given by 

ei = -ibg s , e 2 = -^f 1 . (2.6) 

b 

In this paper, we focus on the b = i case, i.e. the self-dual background e x = — e 2 = g s . The 
momentum conservation condition then reduces to 

2 

m i + m OQ + 2g s N = 0. (2.7) 

i=0 
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This matrix model is expected to reproduce the results of 577(2) gauge theory with Nf = 
4. More precisely, as we will see below, the matrix integral together with the overall factor 
(1 — q) 2 s'i in ( 12. 4 p corresponds to the SU(2) gauge theory. Note that the factor (1 — q) 
is the inverse of the U(l) factor discussed in p]. (See appendix [Al) 

Another point is that the Coulomb moduli parameter a of the gauge theory is identified as 
the filling fraction g s Ni, where Ni is a number of screening operators inserted into the same 
contour in Dotsenko-Fateev representation. For the four-point conformal block we introduce 
N\ and N2. The overall condition (12. 7p reduce these two degree of freedom to one which 
corresponds to the Coulomb modulus of 577 (2) theory. 

The parameters m« are the masses associated with the S77(2) 4 (c 5*0(8)) flavor symmetry. 
These are related to the masses of the hypermultiplets as 

"1l = -(//! + fJL 2 ), ™2 = 2^ 3 + ^ 

1 1 

rn 0O = -(fj, 1 - fJv), m = -(/x 3 - A*4j- (2.8) 

The matrix models associated with gauge theories with Nf = 2,3 are obtained by taking 
the decoupling limit of heavy flavors [32] . By taking a limit of \x± — > oo while keeping \i±q = A3 
fixed, the matrix model action becomes 

W(z) = ^\ogz + m 1 \og(z-l)- 1 ^. (2.9) 

with the following condition: 

mi + moo + /i 3 + 2g s N = 0. (2. 10) 

m l A 3 

The pref actor in front of the matrix integral (12 .4p reduces to e *^ in this limit. This is 



identified with the (inverse of the) U(l) factor of Nf = 3 theory (see appendix IA.2I) 



In order to obtain the Nf = 2 matrix model, we further take the limit \Li — > 00 while 
keeping fi2^3 = ^2 fixed. The dynamical scale of this gauge theory is given by A2. After 
rescaling z — > j&z, the action (I2.9p becomes 

W(z) = fi 3 \ogz-^(z + -) . (2.11) 



2 V z. 

The mass relation reduces in this case to \i\ + ^ + 2g s N = 0. The prefactor becomes simply 



A 
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3 Planar free energy and prepotential 

In this section, we will evaluate the planar free energy of the matrix models introduced above. 
In [12], we have shown that the free energy of these models satisfies the identities known in 
Seiberg-Witten thery [THJ, [19j [20] • Here, we will evaluate the free energies explicitly and compare 
them with the instanton expansions of the prepotentials at lower orders. The computation is a 
bit simpler than in the Seiberg-Witten theory where both the A and B cycle integrals have to 
be computed j2TJ [22] . Here we only have to compute the A integral. 

We first consider the matrix model for SU(2) gauge theory with Nf = 2 in next subsection. 
Then, we will analyze the cases of Nf = 3 and 4 theories in turn. 



The matrix model action corresponding to the SU(2) gauge theory with Nf = 2 is given by 
(12. lip . For simplicity, we will omit the subscript 2 of the dynamical scale A 2 below. There are 
two saddle points determined by the classical equation of motion: 



3.1 577(2) gauge theory with N f = 2 




(3.1) 



These lead to the two-cut spectral curve. 
The planar loop equation reads as usual 




(3.2) 



where the resolvent is defined by 



R{z) = <£ 



z — \j 



>■ 



(3.3) 



l 



The function / is given by 




(3.4) 



Coefficients c\ and c 2 are defined as 




(3.5) 



In the formula for c\ we have used the equations of motion W'(Xj)) = 0. 
Then, the spectral curve x 2 = (2(R(z)) + W'(z)) 2 = W'(z) 2 + f(z) is given by 




(3.6) 
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This is similar to the curve obtained in [23] . The differential one form is identified with A m = xdz 
which has double poles at t = and oo with residues /13 and fi\. Note that the parameter c 2 
corresponds to the variable u in Seiberg-Witten theory. 
We evaluate the filling fraction as 

9sNi = ^-£ A m (c 2 ), (3.7) 

where C\ is a cycle around one of the cuts in the curve. This integral is identified with the 
Coulomb moduli a in the gauge theory and we invert the above relation to solve the unknown 
parameter c 2 . 

Let us compute the free energy of our model defined by 



\I=1^ / KJ 



The starting point is the formula for the A derivative: 

dF m g 



OA 2 



(3.9) 



The expectation value Xj) = (trM) in the second term can be determined by studying the 
large z behavior of the resolvent: R(z) = — \{W'{z) — x) w Ss ^- + 9s ^ r 2 A/ ^ + . . . 

g s (tvM) = -^ K (c 2 -fi 2 1 + ^). (3.10) 

Therefore, we obtain 

A^ = \(2c 2 -^ + ^). (3.11) 

Our remaining task is to determine c 2 in terms of g s Ni by using (13.71) . and this leads to the 
explicit form of the free energy. 

To derive c 2 , let us consider the derivative of (13. 7h with respect to c 2 : 

= / 1 * , (3.12) 

where P4 is the polynomial of degree 4: 

„o 4 . 4/xi 3 4 2 A 2 2 4/i 3 

P 4 ( 2 ) = 2; +^-2; +^2-(/i 3 + c 2 - — )z + — 2 + 1. (3.13) 

It is easy to transform this polynomial so that (I3.12p becomes the standard elliptic integral of 
the first kind. In the following, we set A = /i| + c 2 — 4r and express the result in terms of A. 
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For simplicity, let us consider the equal mass case: [i\ = /X3 = m in the following. In this 
case, by the transformation z = and rescaling of t, the integrand of the right hand side of 
(13.121) can be brought to the standard form 

a/2 dt . 

—= =, 3.14 

^S + (A 2 + 4mA + 2A) ^/(l - t 2 ){l - k 2 t 2 ) 

where k 2 = S-/S+ and 

S± = -^ (-3A 2 + 2A ± AV8A 2 - 1GA + 16m 2 ) . (3.15) 

A 2 + 4mA + 2A\ J v ' 

Then, we can identify the integral ( 13. 12ft in terms of the hypergeometric function: 
. .d{g s N x ) 2V2 f 1 dt 

47TZ- 



dA v /,S + (A 2 + 4mA + 2A) J 1/k y/(l - t 2 )(l - k 2 t 2 ) 

-_F(-,-,l;l-k 2 ). (3.16) 



^/S + (A 2 + 4mA + 2A) 2'2' 



where we have used f, 1 ,, . dt = iK'(k) = iK(k') with k' 2 = 1 — k 2 . We express the 

Jl / k v /(l-t2)(l_fc2 t 2) \ J \ J r 

right hand side as a small A expansion which corresponds to the instanton expansion in gauge 
theory. (Note that k 2 = 1 + O(A).) After integrating over A, we obtain 

r-( m 2 A 2 - 6Am 2 + 15m 4 A , 5m 2 (3A 2 - UAm 2 + 21m 4 ) A , 

2q s N 1 = VA 1 -A A 4 -A 6 

AA 2 64A 4 256A 6 



15 (A 4 - 28A 3 m 2 + 294A 2 m 4 - 924Am 6 + 1001 



m 



16384A 8 

Then, we invert this equation and solve for A: 



A« + 0(A iU ) . (3.17) 



9 m 2 . 9 a 4 — 6m 2 a 2 + 5m 4 , A vr^i'ha 4 — 14m 2 a 2 + 9m 4 ) , K 

A = a 2 H -A 2 H A 4 H - -A 6 

2a 2 32a 6 64a 10 

5a 8 - 252m 2 a 6 + 1638m 4 a 4 - 2860m 6 a 2 + 1469m 8 , „ ,„ /Aln , 

+ r-j A 8 + CA 10 , 3.18 

8192a 14 v ; ' K ' 

where we have introduced a = 2g s N\. Finally, we substitute this into A3. 11 j) and integrate by 

A to obtain 

, 9 9N , a 2 + m 2 , 9 a 4 — 6m 2 a 2 + 5m 4 , , m 2 (5a 4 — 14m 2 a 2 + 9m 4 ) . K 
4F„ = 2 - m") log A + -^A* + ^ A< + ^ ^ 

5a 8 - 252m 2 a 6 + 1638m 4 a 4 - 2860m 6 a 2 + 1469m 8 A „ , nl A in , 

H n A 8 + 0(A 1U ). (3.19) 

32768a 14 K J K ' 

This agrees with the U(2) gauge theory prepotential with a = (a, —a) obtained from the 

Nekrasov partition function (|A.16[) or from the Seiberg-Witten theory [22]. (The first term is 

the one-loop part and the others are the instanton part.) Together with the prefactor e ^ we 
see that the full free energy is the same as that of SU (2) gauge theory. 



6 



3.2 SU(2) gauge theory with N f = 3 

Next, let us consider the matrix model corresponding to the gauge theory with Nf = 3. The 
matrix model action is given by ( 12.9ft . We will omit the subscript 3 of the dynamical scale 
A 3 from now on. As in the previous subsection, there are two saddle points in the classical 
equation of motion. In the planar limit, the loop equation leads to the spectral curve x(z) 2 = 
W'(z) 2 + f(z) where f(z) is written as 

/(*) = - + T^T + * • ( 3 ' 2 °) 
z z — 1 z l 

with coefficients 

— <*(Efe + ^)). <*- ^(E^). * = -%.(e£).<3.21) 

We can easily see that C\ + c 2 = due to the equations of motion. 

The one form defined by A m = x{z)dz has a double pole at z = and a simple pole at 
z = 1 and oo with residues /i 3 , mi and moo, respectively. The residue at z = 00 gives a further 
constraint on q: 

c 2 + c 3 = m 2 00 -{fi 3 + m 1 ) 2 . (3.22) 

This condition together with the relation c\ + c 2 = leaves only one of the parameters inde- 
pendent. Let us choose c 3 to be independent. 

It is then related to the filling fraction by the integral 

Amg s Nx = I A m (c 3 ). (3.23) 

JCi 

For completeness, let us write down here the explicit form of the curve x 2 = iz ^-\) 2 with 

P 4 (z) = Am^z 4 + 4((/i 3 + mi)A + m\ - nl ~ ""4 - c 3 )z 3 

+ (A 2 - 8A/i 3 + 4^3 - 4Am x + 4c 3 )2 2 - 2A(A - 2// 3 )z + A 2 . (3.24) 

It is convenient to introduce the notation B as 

B = c 3 - /i 3 A + /1 2 . (3.25) 

The polynomial is then rewritten as 

P A (z) = 4m^ 4 + 4(Ami + m 2 - - 5)^ 3 + (A 2 - 4A(// 3 + mi) + AB)z 2 

-2A(A-2/i 3 )2 + A 2 . (3.26) 
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Let us consider the free energy of this matrix model. From the definition, its derivative in 
A is written as 

^-f(E^S-4^ A -"3>. P-*) 
In order to determine B we take a derivative of (13.231) with respect to B: 

iJJMA = _ I (3 . 28) 

For simplicity, we consider the case where /X3 = m and m\ = = in what follows. In this 
case, P 4 becomes a polynomial of degree 3: 

P 3 (z) = (2 - l)(-4Bz 2 + (A 2 - 4Am)2 - A 2 ). (3.29) 

After a change of variable (first shifting z — > z — p and then rescaling as z = Qt), we obtain 

P 3 (z) -45Q 2 (1 +p) x t(l - t)(l - fc 2 *), (3.30) 

where 

fc2 = lf^ P = K"i^ (A - 4m) + g )' ^ = ^V / ( A - 4m ) 2 -16P. (3.31) 
As a result, (I3.28P becomes 

A J{9sN l )_ 1 <ft 



dB ^- B (l+p) J y/t(l-t)(l-kH) 

F(H,l;k 2 ). (3.32) 



-B(l+p) v 2'2 

By expanding the hypergeometric function and then integrating over B, we obtain 

/ — / mA 1 9n , 9 m , 9 . ,o 

2^ = ( 1 + — - — 2 (B + 3m 2 )A 2 + ^(5m 2 + B)A* 



(3B 2 + 30m 2 B + 175m 4 ) A 4 -(9B 2 + 70m 2 B + 441m 4 ) A 5 



163845 4 v ' 655365 s 

1 



(5B 3 + 105m 2 B 2 + 735m 4 5 + 4851m 6 )A 6 + C(A 7 ) . (3.33) 



10485765 6 

We invert this equation for B, 

„ 9 mA m 2 + a 2 . 9 a 4 — 6a 2 m 2 + 5m 4 m //lw99 ,^ , 

5 = a 2 + —A 2 + A 4 + - 9a 4 + 70m 2 a 2 + 441m 4 A 5 

2 32a 2 8192a 6 16384a 8 V ; 

2 



+ 7^rr( 1 85a 4 + 1946m 2 a 2 + 15885m 4 )A 6 + C(A 7 ), (3.34) 

262144a iU 
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where we have defined a = 2g s N\. Finally, by substituting this into (j3.27|) . we obtain 

, „ . o 9 . , mA m 2 + a 2 . 9 a 4 — 6m 2 a 2 + 5m 4 , , 
4F m = (a* - m») log A + — + -^A* + ^ A' 

m 9a 4 + 70mV + 441m 4 , «, m 2 185a 4 + 1946m 2 a 2 + 15885m 4 , „ _ /a7 , 

H A 5 H A b + 0(A 7 ). 

2 14 x5 a 8 2 19 x3 a 10 1 J 

(3.35) 

Term with log A is the one-loop contribution. Remaining terms agree precisely with the prepo- 
tential obtained from the Nekrasov partition function (IA.13j) . 



3.3 577(2) gauge theory with N f = 4 

We now consider the matrix model with the original action (12. 5p . The planar loop equation 
R(z) = — | (w'(z) — \JW'{z) 2 + J[z)j involves a function f(z) which now has a form f(z) = 
Y^,=q Parameters {q} are given by 



c 



-4o s m (5^T-), ci = -Ag s nii{ 1 ), c 2 = -Ag s m 2 { V" t- - — ). (3.36) 

j A/ ^y" A/ — 1 ^ A/ — g 



By studying the behavior of loop equation at large z we find that the parameters obey 



2 



^Ci = 0, ci + gc 2 = m 2 ^ - (J^m^ 2 . (3.37) 

i=0 i=0 

By eliminating ci and c 2 , the spectral curve becomes 

2 _ ^(z) , , 

x - z 2 (z-i) 2 ( Z - q r 1 j 

where P4 is the following polynomial of degree 4 

P 4 (z) =m 2 XD z 4 + ( - (1 + g)(m 2 x> + mo) + (1 - - m 2 ) - 2m (gm 1 + m 2 ) + gc jz 3 

+ (<? m L + (1 + 3g + g 2 )mg + (1 - q)(m\ - qm\) + 2(1 + q)m (qm 1 + m 2 ) 
-(1 + q)qc J z 2 + ^ - 2g(l + g)mp - 2q 2 m mi - 2qm m 2 + g 2 c j 2; + q 2 m 2 Q . (3.39) 

The meromorphic one form xdz has simple poles at z = 0,1, g and z = 00 with residues 
mo,mi,m 2 and moo. 

Again, we consider the derivative of the free energy: 

3p / 1 \ 

9 s m 2 {tT- -) = m2 R( z )\ z=q . (3.40) 



dq \ q — M 
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This can be easily computed by expanding the resolvent at z — q, R(z) = + 0(z — q). 
Then, we obtain a simple expression for the free energy 

^ = 7=4(1^) (<|>>*-"4- *)■ ("D 

In the last equality we used the relation (I3.37p . 

In what follows, we consider the simple case where all the hypermultiplet masses are equal 
to m: i.e. mo = = and mi = m2 = m. In this case, the polynomial is reduced to degree 
3: P?X Z ) = Cz{z — z + )(z — Z-), where we have introduced C = c Q q and 



z± = ±\l + q-(l- qf^ ± (1 - q))jl - 2(1 + q)^ + (1 - q)^j . (3.42) 

By taking the C derivative of xdz, the holomorphic one form becomes 

d 1 dz 2 z 2 _ 

^—-xdz = - . — , k = — . (3.43) 

dC 2^/Uz^ y/z(l - z){l- k 2 z) q V ; 



The remaining calculation is similar to those considered in the previous subsections. That is, 

o 

m 
C 



2 

we first evaluate the period integral of the above one form. Then by expanding in and 



integrating over C, we obtain 



2ig.N 1 = VC [ h (q) - h^q)— - — - + O(-) ) , 



,2 



where hi(q) are the expansion coefficients of the period integral in and depend only on q. 
ho(q) is for the theory with massless flavors: 

h (q) = l + l - q+ ^ + + iffV + 0(q 5 ). (3.45) 

ow V 6V 256 y 16384 ^ w; v ; 

Lower order expansions of hi, hi and h^, are given by 

^ = \ + \ q+ w/ + 5iV + ^/ + °^ 

h ^ =l + ^ q+ §2 q2 + ms q3 + isi^^ 4 + ° {q% 

, . . 5 125 1125 2 125 3 125 4 5 . 

h 3 (q) = 1 OH q 2 H o H o 4 + C(g 5 ). (3.46) 

3W 16 64 y 1024 y 4096^ 262144 y w 7 v ; 

Solving for C, we obtain 

1 2h 1 (q)m 2 2h (q) h 2 (q) - 3/ii(g) 2 m 4 



C = a 2 



i 



h (q) 2 h (q) a 2 3 a 

+ lQMg)hi(g) 3 - iQfeo(g) 2 fei(g)fe 2 (g) + 2/i (g) 3 /i 3 (g) m 6 + (3 47) 



10 



where a = 2ig s N\. By substituting the above expression into (13.411) and integrating over q, we 
finally obtain the Nf = 4 free energy 

. o 9N , a 4 + 6a 2 m 2 + m 4 13a 8 + 100m 2 a 6 + 22m 4 a 4 - 12m 6 a 2 + 5m 8 „ 
AF m = (a 2 - m 2 ) log q + — q + 



2a 2 * 64a 6 

23a 12 + 204m 2 a 10 + 51m 4 a 8 - 48m 6 a 6 + 45m 8 a 4 - 28m 10 a 2 + 9m 12 



g 3 + 0(g 4 ).(3.48) 



192a 10 

This perfectly agrees with the instanton partition function flA.9j) . 

Finally, we make a brief comment on the massless theory. In this case, the expression for 
C simplifies and becomes C = a 2 /h^{q) where h (q) is (I3.45p . Thus, it is easy to derive 

AF m = a 2 ( log q - log 16 + -q + —q 2 + — q 3 + -^q 4 + ^^q 5 + 0(q e )] , (3.49) 
V 2 64 192 32768^ 81920^ W ') ' V ; 

where we have added the one-loop contribution —a 2 log 16. Note that this can be written as 
4F m = a 2 log qiR where q and qiR = e 2ntTIR are related by 

1 = M^JI = 16( ^ R - 128 ^R + 704 ^ " 3072 ^ + 11488 ^R + • • • • ( 3 - 5 °) 

as already discussed in [251 12S1 [H |271 [71 [12]. Thus the theory appears classical in terms of IR 
coupling constant r IR . 



4 Matrix model and Quiver gauge theories 

In this section, we study matrix models which describe M = 2 SU(2) quiver gauge theories. 
First of all, we consider a matrix model describing SU(2) linear quiver gauge theory where 
each gauge group has a vanishing beta function [28J. Then by taking its decoupling limit, we 
propose models for asymptotically free gauge theories in subsection 14.11 

According to the AGT conjecture, SU(2) n ~ 3 linear quiver gauge theory is related to the 
n-point conformal block of the Liouville theory, which is represented by the trivalent graph [29] 
as in Fig [TJ As seen in section [21 the Dotsenko-Fateev representation of the conformal block 
suggests a matrix model with the following action [TT] : 

n-2 

W{M)=^2milog{M-ti), (4.1) 

i=0 

where t — and t\ = 1. Other parameters tj = Ylk=i Ik (i — 2, . . . ,n — 2) describe complex 
structure of the n-punctured sphere. Note that we also have the prefactor as in (12. 4p 

From the gauge theory perspective, the parameters qt are related to the gauge coupling 
constants qu = e 2niTk of the gauge group SU(2) n ~ 3 . For n — 4, this reduces to the matrix 



11 



rn\ J712 m 3 m n-2 



H 12 



moo tuq 



Figure 1: 



model which we studied in subsection 13.31 Parameters m and m n _ 2 are related to the mass 
parameters of two hypermultiplets in the fundamental representation of the SU(2) at one end 
of the quiver: m n _ 2 = |(/i3 + fi^) and m = |(/i3 — fii). Also, the masses of the bifundamental 
hypermultiplets are identified with m« (i = 2, . . . , n — 3). Finally, the masses of two fundamental 
hypermultiplets of the other end of the quiver are related to m\ and as m\ = |(/ii + /i 2 ) 
and moo = §(^1 — ^2)- The mass parameter is introduced by the following condition: 

n-2 

m + moo + 2g s N = 0. (4.2) 
The critical points are determined by the equation of motion 

EsTV^Ej^Aj^ (43) 

i=0 J 1 J(^/) J J 

If we ignore the second term, we obtain n — 2 critical points e p (p — 1, 2, . . . , n — 2). Let each 
iV p (p = 1, 2, . . . , n — 2) be the number of the matrix eigenvalues which are at the critical point 

Cp. 

We take the large iV limit with mass parameters {mj} and filling fractions \v v = g s N p } 
being kept fixed. Since this is one matrix model, the loop equation is still the same as in the 
previous cases (13. 2p 



^^^'l-^Lg-S-E fW . (4.4) 

We note that a polynomial Z(t) is of degree n — 3, since the leading term vanishes due to 
equations of motion. 

Finally, we define the meromorphic one form A = x(z)dz as 

x(z) 2 = {2{R{z)) + W\z)f = W\zf + f(z). (4.5) 
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4.1 Matrix model for asymptotically free quiver gauge theory 

The matrix model corresponding to asymptotically free quiver gauge theory can be obtained 
by taking the decoupling limit as in section [2] Only possible limits which does not spoil the 
condition f!4.2[) is the case where /i 2 (= m i ~ m oo) or rn n _ 2 — m ) is taken to infinity. 
For the sake of illustration, let us consider the n = 5 case with the action 

3 

W(z) = J2 m i^g(z-t i ) J (4.6) 

i=0 

where t 2 = q\ and t 3 = q\q 2 . This corresponds to S77(2)i x SU{2) 2 quiver gauge theory whose 
gauge coupling constants are qi and q 2 . We first take a limit fi^ — > oo with /x 4 g 2 = A fixed. In 
this limit, we obtain 

W{z) -> /i 3 log z + m i - **) - y^- ( 4 - 7 ) 

i=l,2 

It is natural to anticipate that this matrix model corresponds to the quiver theory of one 
fundamental matter coupled to the second gauge group SU{2) 2 and two fundamental multiplets 
are coupled to the first gauge group. The relation of the mass parameters H4.2f) becomes 
/i3 + X^=i,2 mi + m °° + %9sN = in this limit. 

By further taking the limit fi 2 — > oo with = A fixed, we obtain from (" 14. T[) 



W(z) ->■ /i 3 log z + m 2 log(z - 1) - ^ - y> ( 4 - 8 ) 

where we have also rescaled z — > q\z. The relation of the mass parameters (14.21) becomes 

/is + m 2 + [i x + 2# s iV = 0. (4.9) 

This matrix model is expected to describe SU(2)i x SU{2) 2 quiver gauge theory with each 
gauge factor coupled to one hypermultiplet. Both of the gauge factors have nonvanishing beta 
functions and the theory is asymptotically free. 

It is possible to generalize this construction to the case with n > 5. A decoupling limit of 
a hypermultiplet at the last end of the quiver is /A4 — > 00 with \i/±q n -z = A fixed. Also, another 
decoupling limit of a hypermultiplet at the first end of the quiver is [i 2 — > 00 with /i 2 gi = A 
fixed. By taking these limits, we finally obtain 

n-3 / i-1 \ ( n ~ A \ 

W{z) = fi 3 log z + m 2 \og{z - 1) + ^Vilog I z- ]Jq k I - - — I JJ q k I , (4.10) 

i=3 \ k=2 I \k=2 J 



with the following relation for the mass parameters: 

n-3 

/i 3 + J2m i + v 1 + 2g s N = 0. (4.11) 

i=2 
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5 Conclusion and discussion 



In this paper we have studied the matrix model proposed to explain the AGT relation and 
interpolate the Liouville and M = 2 SU (2) gauge theories. We have explicitly evaluated the 
free energy of the matrix models describing SU(2) gauge theory with Nf = 2,3,4 flavors and 
have shown that they in fact reproduce the amplitudes of Seiberg-Witten theory. Our analysis 
is limited to the large N limit and it is very important to see if our results can be generalized 
and reproduce full Nekrasov partition functions. There is already an interesting work in this 
direction [301 EI] and we hope that we can report further results in future publications. 
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Appendix 

A Nekrasov partition function 

The instanton partition function of J\f = 2 U(2) gauge theory with Nf = 4 is expressed as a 
sum over all possible Young tableaus parametrized as Y = (Ai > A2 > . . .) where A^ is the 
height of the £-th column [2U [1] : 

^inst — ?' y ' ^vector ( Q ) ^)Zantifund(<J, Y, /il)Zantifund(o, Y, jl-z) Zf un d (a, Y, — /X3) Zf un d (a, Y, —^4). 

(A.l) 
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Here 



Z vcctor (a, Y) = Yl II ( a v ~ 6 ^ L yA s ) + ^(A Yl (s) + 1)) 
ij=i,2 sgy, 

II (%' + £ i L ^ (*) - e 2(A Fi (t) + 1) + e + ) _1 , 



x 



^fund(a, = ] [ (a; + ei(^ - 1) + e 2 (m - 1) - fi + e+) , 

i=i,2 seVi 

^antifund(a, F, /i) = JJ JJ (a, + e x {l - 1) + e 2 (m - 1) + //) , (A.2) 

1=1,2 sGY 

and e + = ei + e 2 and For a box s at the coordinate (£,m), the leg-length 

-ky(s) = X' m — i and the arm-length A Y {s) = \& — m where \' m is the length of the m-th row. 
The minus signs of the masses in Zf un d are due to the convention. 

In order to derive the expression for 577(2) gauge theory, we set the Coulomb moduli as 
a = (a, —a) which gives 

■^vector (fli 

Y)=] [ (2ac% - e 1 Ly j ( 8 ) + e 2 {A Yi {s) + 1)) 1 

i=l,2 s&i 

x \[ (-2a5ij + e 1 L Yj (t) - e 2 (A Yi (t) + 1) + e+) _1 , 
Zfund(a, Y,fj)= ] [ (a6i + ei(£-l) + e 2 (m - 1) - fi + e+) , 

i=l,2 s&Yi 

^antifund(a, ?,//)=] [ (a^ + <i\{£ ~ 1) + 6 2 (m - 1) + /i) , (A.3) 

i=l,2 seFi 

where we define 5i = +1 and <5 2 = —1, and 

{0 for % — j, 
1 for i = 1 and j = 2, (A.4) 
— 1 for z = 2 and j = 1 . 

Then, the SU(2) and 17(2) partition functions are related by the U(l) factor as pointed out in 

a- 

Note that this expression differs by a minus sign in front of (/X3 + 114) from the one of [I]. As 
argued in pQ, the SU(2) partition function is invariant under "flips". These flips are reduced 
in the self-dual case t\ = — e 2 = h to 

a-^-a, fix ±fl 2 -»■ -(fix ±/i 2 ), yU 3 ±yU 4 ^ -(/i 3 ±yU 4 ). (A. 6) 
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Gauge theory prepotential can be obtained in the limit where the deformation parameters 
go to zero (with a fixed ratio £1/62): 

Jlnst = lim (-eie 2 )logZ inst . (A.7) 

In the self-dual case, 577(2) gauge theory prepotential is written as 

TZ i2) = ^(logZ inst -logf u U) 

= + 2(^1 + ^2)^3 + /x 4 )log(l - g). (A.8) 

To compare with the free energy of the matrix model, we present an expansion of J-"i ns t for the 
equal mass case /Xj = m 

a 4 + 6m 2 a 2 + m 4 13a 8 + 100m 2 a 6 + 22m 4 a 4 - 12m 6 a 2 + 5m 8 „ . , 

Jw= ^ ?+ sis* « (A ' 9) 

23a 12 + 204m 2 a 10 + 51m 4 a 8 - 48m 6 a 6 + 45m 8 a 4 - 28m 10 a 2 + 9m 12 , ml 

+ 192^ " + 0(9 >■ 

A.l U(2) gauge theory with Nf = 3 

Let us consider Nekrasov partition function of the theory with Nf = 3. This can be obtained 
from the above partition function by taking a limit \l± — > 00 with /i 4 g = A 3 fixed. In the 
/c-instanton part the only factor which contains /X4 is 

Z hmd (a, Y, -114) = II II ( a6i + - 1) + e 2 (m - 1) + /i 4 + e+) . (A.10) 

i=l,2 s€Yi 

When combined with fc-instanton factor q k , this gives the leading contribution A3 and the other 
contributions are suppressed in the limit. Therefore, we obtain 

^inst = M ^vcctor( a ; ^)-^antifund(°; ^ i / i l)^antifund( ( 2 ) Y, fl 2 ) Z{ und (a,Y , — . 

c*i,y a ) 

(A.ll) 

The t/(l) factor reduces to 

f U(l) f U(l),N f =3 = exp /_ (^l+/i 2 )A 3 ^ (A 12) 

V 2eie 2 / 

In the simple case of /13 = m and /ii = /i 2 = which we considered in subsection 13.21 the 
prepotential of the gauge theory is given by 

_jVf =3 1 a 2 + m 2 9 a 4 — 6m 2 a 2 + 5m 4 . , 

= + + ^ A 3 + ^( A 3)- (A.13) 
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A. 2 17(2) gauge theory with N f = 2 

We can further take a limit where /12 — > 00 while keeping //2A3 = A 2 , fixed. In this limit, the 
partition function becomes: 

^inst = y '^vcctor(o, Y) Z antifund (d, Y, A t l)^fund( ) ~~ A^)) (A. 14) 

and the U(l) factor is reduced to f u ^ — > exp ^— SU(2) prepotential is given by 

T SU(2),N f =2 _ N f =2 A| , A 1 s 

inst — S inst {A. 10) 

For the equal mass case with [x\ = [13 = m, lower terms of instanton expansion are given by 



^.N f =2 a 2 + m 2 t9 a 4 - 6a 2 m 2 + 5m 4 . . m 2 (5a 4 - 14a 2 m 2 + 9m 4 ) , fi ^, iSn 
^ =^^ A2 + 64^ A2 + ^ 192^ ^AS + OCA-). 



(A.16) 
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